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A possibility of a metal-insulator transition in molecular conductors has been studied for 
systems composed of donor molecules and fully ionized anions with an incommensurate ratio 
close to 2:1 based on a one-dimensional extended Hubbard model, where the donor carriers 
are slightly deviated from quarter filling and under an incommensurate periodic potential from 
the anions. By use of the renormalization group method, interplay between commensurability 
energy on the donor lattice and that from the anion potential has been studied and it has been 
found that an "incommensurate Mott insulator" can be generated. This theoretical finding will 
explain the metal-insulator transition observed in (MDT-TS)(Aul2)o.44i- 
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Most of the conducting molecular crystals are realized 
by combining two kinds of molecules, A and B, with 
commensurate composition ratios. Typical examples are 
the well-studied 2:1 compounds, i.e., A^B, which show 
a variety of phases such as Mott insulator, charge order, 
superconductivity and so on. 1 ' 2 In the compounds, the 
molecule B is usually fully ionized either as —1 or +1 
to form a closed shell, and as a consequence the energy 
band formed by HOMO or LUMO of A is quarter-filled 
as a whole in terms of holes or electrons, respectively. 

Recently, molecular conductors with incommensurate 
(IC) composition ratios close to 2:1 have been synthe- 
sized based on new donor molecules. 3-5 (MDT-TSF)A„ 
and (MDT-ST)A n (X = I 3 , Aul 2 or IBr 2 , n ~ 0.42 - 
0.45) show metallic behavior and undergo superconduct- 
ing transition at about T c = 4 K at ambient pressure. 3. 
In contrast, in (MDT-TS)(Aul2)o.44i a metal-insulator 
(MI) crossover occurs where the temperature dependence 
of resistivity displays a minimum at T p = 85 K. 5 In ad- 
dition, an antiferromagnetic transition takes place at Tn 
= 50 K. When pressure is applied to this compound, T p 
decreases and the superconducting phase appears above 
P c = 10.5 kbar (T c = 3 K). All these compounds are 
isostructural with alternating donor and anion layers. 
Since the anions are fully ionized as A" - as in the 2:1 
compounds, 3 the electronic properties can be attributed 
to the donors with the IC band-filling slightly larger than 
3/4 for the HOMO bands. The extended Hiickel scheme 
predicts two-dimensional (2D) Fermi surfaces which are 
similar to each other. 4, 5 It should be noticed that the 
anions in these compounds are not randomly distributed 
in the layers, but are found by the X-ray scattering ex- 
periments to form regular IC lattices with a different 
periodicity from the donors. 3,4 

The metallic state observed in these compounds is nat- 
urally expected from the IC band-filling since the sys- 
tem would avoid insulating states due to strong cor- 
relation such as Mott insulator or charge order. On 



the other hand, it is difficult to understand the strong- 
coupling nature of the insulating ground state in (MDT- 
TS)(Aul2)o.44i, deduced from T p ^ Tn, which is to be 
explored in this Letter; if it is the weak-coupling spin- 
density-wave state due to the nesting of the Fermi sur- 
face, T p — Tn would be expected. 1 

We consider a one-dimensional (ID) model in order to 
capture the essence of the MI transition in a more con- 
trolled way than considering a 2D model relevant to ex- 
periments. Our ID model consists of Nr, donor molecules 
coupled with N anions both forming regular lattices, as 
shown in the inset of Fig. 1. The donors are modeled by 
the ID extended Hubbard model, known to be relevant 
for typical 2:1 systems, 6,7 and the small potential from 
the anions is added, which is crucial for the insulating 
state to appear. The Hamiltonian is written as follows, 
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where t, U and V are respectively the transfer energy 
between the nearest-neighbor donor sites, the on-site re- 
pulsive interaction and the nearest-neighbor repulsion; 
the creation operator at the j-th site with spin s=± is 
denoted as c\j s , n,j, s = ct g Cj tS and rij — Y2 s n J> s - Since 
the fully ionized anions form a regular lattice, the an- 
ion potential at the j-th site, Vj, can be expressed as 
v 3 = N Z 1 Em=-oo v(mQ)e imQja where Q = 4fc F , k F = 
Tin/ (2a) is the Fermi wavenumber, n = N/Nl is the car- 
rier density (we take the hole picture in the following) in 
the donor chain and a is the spacing between donor sites. 
In the following, we consider only the relevant ±Q com- 
ponent of the potential, v(±Q) = vne ±lx . This can lead 
to a gap, 2i>o, at ±2/cp in the non-interacting band, which 
we assume to be small compared to the band width. Then 
the system becomes effectively half-filled in reference to 
the IC lattices as is seen in Fig.l. 
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Fig. 1. Energy dispersion in the presence of the anion potential 
i>0 = O.lt where the occupied onc-particle states are expressed by 
the thick curve. The figure is written in the case of n = 0.33 to 
clarify the characteristics of the present model. Inset: a schematic 
representation of our model. 



We derive an effective Hamiltonian for low energy scale 
in terms of phase variables following the just quarter- 
filling case. 7 To the lowest order of the normalized an- 
ion potential S = v / {e(3kp) — e(kp)} with e(K) = 
— 2t cos Ka, the phase Hamiltonian is obtained as H c g = 



TLp+TL a +TL 



pi j. Here Tip, Ti a and TLp a are respectively the 



charge part, the spin part and the term mixing both de- 
grees of freedom. The spin part, 7i a has the same form 
as that of the ID Hubbard model, so the spin excita- 
tion becomes gapless. 8 The term H pa is expressed by 
the product of the non- linear terms seen in H p and H a , 
and then has a larger scaling dimension. Hence we can 
neglect it. 9 Therefore the properties of the charge degree 
of freedom are determined by H p , expressed as 
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Here qo — 2ir/a — 8fcp = 2n(l — 2n)/a is the misfit param- 
eter, v = 2tasinfcpa, v p = v^B p A p and K p = yj B p /A p 
with A p = 1 + (,g 4 || + g 4± + g 2 \\ + g 2 ± - 9i\\)/{nv) and 

52± + 9i\\) / {nv) . or 1 is the 
a). The interaction parameters 
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ultra-violet cut-off (a 
are written as 
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where v(3kp) ~ 2tasin3fc F a. 

In eq.(2), there are three non-linear terms. First, the 
half-filling Umklapp term, , is generated by the anion 
potential 5 because the band is effectively half-filled as 
seen in Fig.l. This can lead to a Mott insulator, as we 
will later show explicitly. We call the state as IC Mott 
insulator since it has a periodicity not matching with 
the donors but with the anions. However, it is not trivial 
whether this IC Mott insulator can be realized, and if 
so, in which condition it is stabilized, in contrast to the 
half-filled Hubbard model where infinitesimal on-site re- 
pulsion stabilizes the Mott insulator. 8 This is because of 
the presence of the other two non-linear terms in eq.(2), 
the "quarter-filling" Umklapp term, 51/4, with the misfit 
which is present even without the anions owing to the 
proximity to a quarter filling on one hand, and the g'^ ± 
term, the combination of both commensurabilities of the 
donors and the anions on the other hand. 

In the present calculation, it is crucial to fix the car- 
rier number at the value determined by the anion density. 
However, it is to be noted that if one evaluates the quan- 
tity based on eq. (2) , it results in a deviation of the car- 
rier number from the value in the non-interacting case, 
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where G' 3 \ = S r 3J_/( 7r ^), G 1/4 = g 1/4 /(2wv p ), and J n (x) 
is the the Bessel function of the first kind. The origin of 
the deviation is the existence of the misfit parameter. 10 
Therefore, we must add the term, —((i/ir) J dxd x 9 p and 
keep AN e to zero. To the lowest order of G' 3 ' ± and G1/4, 
the chemical potential fj, is given as, 



q„a = -4K p G'£ 



-8K n G 2 



J a a a 
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where q p = 4K p fi/v p . 

To determine the low energy behavior of this effective 
Hamiltonian, we derive the renormalization group (RG) 
equations by rewriting the action S p corresponding to 
the Hamiltonian, H p — (n/n) J dxd x 9 Pl as 

; f d 2 r{(dj p ) 2 + (d y 9 p ) 2 } 
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where 9 P = 9 p -q p x/2, G 3 ± = g 3 ±/{irv p ), q 1/4 = q Q -2q fl 
and q' 3 ' = <7o — q^- The condition, AN e = 0, leads to the 
following self-consistent equation, 
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where Ji(x) = sgn(x)Ji(\x\). Eqs.(16) and (17) lead to 
the following RG equations, 
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Fig. 2. The solutions of the RG equations, K p and (91/4 — q'l)ct ( 
the misfit parameter of the G3T-term ) for U /t = 3.0, V/t = 0.0 
and n = 0.441. The cases of <5 = 0.001 and 0.005 are denoted by 
the solid and dotted curves, respectively. 
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Eqs.(18)-(23) are obtained for the condition of the ac- 
tion, eq.(16), being invariant under RG transformation, 
whereas eq.(24) is obtained from the condition of the 
chemical potential, eq.(17). From eqs.(19), (20) and (24), 
it is shown that the quantity q , whose dimension is 
(length) -1 , is scaled as q (l) — q a e l . This shows the fact 
that the carrier number is indeed conserved without any 
effects from the interaction. 

Typical flows of the RG equations are shown in Fig. 2 
for U/t = 3.0 and V/t = 0.0 where the carrier num- 
ber is fixed as n = 0.441 taken from the actual mate- 
rial (MDT-TS)(AuI 2 ) .44i- In the case of S = 0.005, K p 
tends to zero implying that the ground state is an insula- 
tor, due to the commensurability in the half- filling, G 3 ±. 
This can be seen in the RG equations since the misfit pa- 
rameter in the G3^-term, (qi/ 4 — q 3 )a = —q p a vanishes 
(see Fig. 2) and then G3j_ affects the renormalization of 
K p through eq.(18) while those in the G1/4- and G 3± - 
terms, 91/40: and q 3 a, tend to 00 and then these effects 
become negligible due to the oscillating behavior of the 
Bessel function. Hence the origin of this insulating state 
is nothing but the commensurate potential of the effec- 
tive half-filling generated by the anion potential. Namely, 
the insulating state is indeed the IC Mott insulator. 

On the other hand, in the case of smaller potential 
due to the anions, S = 0.001, a metallic state with fi- 
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Fig. 3. The phase diagram on the plane of <5 and U/t in the case 
of V/t = and n = 0.441. 



to 




Fig. 4. The phase diagram on the plane of 8 and V/t in the case 
of U/t = 3.0 and n = 0.441. 

nite value of K p is realized. Here in eq.(18) the effects 
of the (j?3_i_-, G% ± - and Gi/4-terms on K p disappear at 
the low energy since all the misfit parameters are diver- 
gent. This metallic state is not realized if we set G1/4 to 
zero. Therefore we can state that the origin of the MI 
transition is the interplay between the different kinds of 
commensurabilities. 

Next, we show ground state phase diagrams as a func- 
tion of the model parameters. First, the phase diagram 
on the plane of U/t and S in the case of V/t = is 
shown in Fig. 3. Since the quantity 5 is proportional to 
Vo/t, the transition from the metallic state to the IC 
Mott insulator occurs when the potential from the an- 
ions increases and/or the band width decreases. When 
U —►00, the present system can be mapped onto a 
non-interacting spinless Fermion system with the Fermi 
wavenumber doubled, as 2fc F . n In this case the insulat- 
ing state is realized by an infinitesimal <5 because a gap 
opens at ±2fcp (see Fig.l), consistent with Fig. 3. 

The role of the Gi/4-term on the MI transition be- 
comes clearer when V is varied. It is because the cou- 
pling constant G1/4 changes its sign when V increases. 7 
We show the phase diagram on the V/t - 5 plane in Fig. 4 
for U/t = 3.0. At V = V c = 0.838t where G 1/4 = 0, the 
IC Mott insulating state is realized for infinitesimal S. 
For V > V c , the absolute value of G1/4 increases again 
and results in a finite metallic region. Therefore, a re- 
entrant transition, metal — > IC Mott insulator — > metal, 
occurs when V increases. Note that there is no quali- 
tative difference between the metallic states in the two 
distinct regions. 

Finally let us discuss the relevance of our results to 
the experiments. The difference of the ground state in 
metallic MDT-TSF and MDT-ST compounds and that 
in the MDT-TS compound undergoing MI crossover can 
be naturally understood as follows. The extended Hiickcl 



scheme provides transfer integrals, i.e., the band width, 
of the MDT-TSF families larger than that of the MDT- 
TS compound, 4 ' 5 which is consistent with our results 
that the decrease of the bandwidth lead to an MI tran- 
sition, as seen in Fig. 3. In our ID model the spin de- 
gree of freedom is essentially that of the ID Heisenbcrg 
model showing no magnetic order. However, in the IC 
Mott insulating state in the actual 2D material we gen- 
erally expect that antiferromagnetic order appears at low 
temperature due to the three dimensionality, as in fact 
observed. 5 In this case, the magnetic ordered moment 
should be large, compared to, e.g., that of the spin- 
density-wave state due to the nesting of the Fermi sur- 
face. 

In conclusion, we investigated the electronic state of 
the one-dimensional extended Hubbard model close to 
quarter-filling under an incommensurate anion potential. 
We found that a transition between the metallic state 
and an incommensurate Mott insulator can occur, whose 
origin is the interplay between the commensurability en- 
ergy generated by the anion potential and that in the 
donor lattice. To the authors' best knowledge this is the 
first theoretical study of a "Mott transition" generated 
by such interplay between different commensurabilities. 
It would be interesting to investigate the critical proper- 
ties of this transition in the actual compounds and com- 
pare with the "usual" Mott transition seen in the typ- 
ical 2D A 2 B molecular conductors, k-(BEDT-TTF) 2 X, 
which is recently attracting interests. 12 
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